Accepted for Publication in ApJ 

Preprint typeset using LAT^^ style cmulateapj v. 14/09/00 



(N 
O 

o 

(N 
oo 



MASS SEGREGATION IN GLOBULAR CLUSTERS 
J. M. Fregeau\ K. J. JosHi^, S. F. Portegies Zwart''"', F. A. Rasio^ 

Accepted for Publication in ApJ 

ABSTRACT 

We present the results of a new study of mass segregation in two-component star clusters, based on 
a large number of numerical A^-body simulations using our recently developed dynamical Monte Carlo 
code. Specifically, we follow the dynamical evolution of clusters containing stars with individual masses 
mi as well as a tracer population of objects with individual masses m2. We consider both light tracers 
{fi = 7712/7711 < 1) and heavy tracers (/x > 1), and a variety of King model initial conditions. In all of our 
simulations we use a realistically large number of stars for globular clusters, N — 10^, but we ignore the 
effects of binaries and stellar evolution. For heavy tracers, which could represent stellar remnants such 
as neutron stars or black holes in a globular cluster, we characterize in a variety of ways the tendency for 
these objects to concentrate in or near the cluster core. In agreement with simple theoretical arguments, 
we find that the characteristic time for this mass segregation process varies as For models with 

very light tracers (/i < 10^^), which could represent free-floating planets or brown dwarfs, we find the 
expected depletion of light objects in the cluster core, but also sometimes a significant enhancement in 
the halo. That is, for some initial conditions, the number density of light objects in the cluster halo 
I increases over time, in spite of the higher overall evaporation rate of lighter objects through the tidal 

^ ■ boundary. Using these results along with a simplified initial mass function, we estimate the optical depth 

f> — ' to gravitational microlensing by planetary mass objects or brown dwarfs in typical globular clusters. For 

l/^ ! some initial conditions, the optical depth in the halo due to very low-mass objects could be much greater 

' than that of luminous stars. If we apply our results to M22, using the recent null detection of Sahu, 

Anderson, & King (2001), we find an upper limit of ~ 25% at the 63% confidence level for the current 
mass fraction of M22 in the form of very low-mass objects. 

Subject headings: clusters: globular — celestial mechanics, stellar dynamics — Monte Carlo: dynamical 
evolution — stars: low mass, brown dwarfs, planetary systems 

I ' 1. INTRODUCTION dynamical evolution of these clusters, including more mas- 

O u 1 1 i ii ii i i • u- i -ii sive main-sequence stars and binaries, as well as primordial 

V-j Globular clusters are thought to contam objects with m 7,^ i,i i i, i t-o t- ■ ry \ p a/t ^/r-ii onnn\ 

» ■ .J r 1.1 1 .1, . 10 Mp, black holes (Portegies Zwart & McMillan 2000). 

ry-i , a very wide range of masses, although even the present- ^ ^ i i i ■ , i ^ ^ <• i . 
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' r VrTv/rnx • i j. ■ ^ -i j j i tcr initially and, as we will show m this paper, signihcant 

• I— I , function (IMl* ) IS even less constrained, as it depends also , , , , , ■ X . 

k>( ii J i -1 r J.1 11 1 J. J -1 1 i- numbers could also have been retained to the present. 

' on the details oi the overall cluster dynamical evolution ^ i i • , ,^ i i 

'^(2 ' nr • ■ p Tj • inn'7\ t j j -(• u- 1,1 -j Some lower mass obiects '--^ 10 ^ — 10 Mr have been 

H (Vesperim & Heggie 1997). Indeed, even tor highly idc- • ; -n- 

' \. 1 , f 1 ■ . , r, detected m globular clusters as companions to milhsec- 

. . . 1 ahzed systems ot unevolvmg point masses, a wide mass , ,. ° /t^, » . , „^r;-, , i 

J- 1,- J -iu 4-u a- J- f J- u J 1 ond radio pulsars (D Amico et al. 2001; lord et al. 2000; 

spectrum, combined with the enects ot two-body relax- „ . , \ ^•i,-, , ; , i rmrr^ , 

^- 1 \ . ■ ^ r ij -11, Freire et al. 2001). Gilhland et al. (2001) used HST to 

ation, leads to a variety ot complex dynamical phenomena , „ , . . . , ,. 

, , -n 1 J 4- J j-i, 4.- 11 /tr • J- search tor transits by giant planets m short-period orbits 
that are still poorly understood theoretically (Heggie et /a-, . r „\ i • , • ;i ; i 
1 inno\ -u f 4- 4-u i, u I- i ( hot Jupiters ) around mam-sequence stars m the Central 
al. 1998). We will refer to these phenomena collectively ^ . r .-m mi .1 ;• 1 
1 a " u J- i j-i, J- j-i, • 1 region ot 47 iuc. 1 hey reported a negative result and con- 
here as mass segregation, but note that they involve a , , . r -.Sm .1 .1 
number of difi^erent processes such as as mass loss through '^I'^^ed that the planet frequency m 47 Tuc must be at least 
ii, i-j 1 u J • J J. an order ot magnitude below that tor the solar neighbor- 
the tidal boundary, energy equipartition and mass stratifi- , , -,t • ;i i • 1 1 •; ,1 • n.i- 
,. 4-u 1 J. i- J • hood. However, m the high-density central region 01 this 
cation, gravothermal contraction, and even, m some cases, , , , , , 1,11 

,1 1 • , 1 -Ti- / -i. Tnan\ clustcr, whcrc the search was conducted, planetary sys- 

gravothermal instabilities (see, e.g., Spitzer 1987). , ,-, 1 , 1 i- ; 1 1 ; -.i 

fj- ii u uij- ij - lui tems are likely to be disrupted by encounters with other 

ihe masses of directly observable stars today m globular , /t^ • r.- 1 1 

1 , w 1 c 1 r, stars and binaries (Davies & Sigurdsson 2001), and the 

clusters cover a relatively narrow range from ~ 1 — 2 Af© uenc of "free floatin " lanets is not constrained b 

(e.g., primordial binaries, blue stragglers, neutron stars; requency o lee oa mg p ane s is no cons rame ^ y 

o -1 inncN J 4- r>? i\^ ^ ^x, c ■ i- A these observations. In contrast, gravitational microlensmg 
see, e.g., Bailyn 1995) down to '--^ 0.1 M(7) at the lamt end , , , , ; , , , , n ;• „ 1 - 

n , 1 • / • 4. 1 1nno^ 1, can be used to try to detect such tree-fioatmg planets 

of the mam sequence (e.g., Marconi et al. 1998). Much b i- 

more massive stars were certainly present earlier in the 
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(Paczynski 1994). Analysis of microlcnsing events toward 
the galactic bulge suggests the presence of a substantial 
amount of lower mass objects in globular clusters (Jetzer, 
Strassle, & Wandeler 1998). More recently, in a pioneering 
study, Sahu et al. (2001) monitored about 83,000 Galac- 
tic bulge stars for microlensing by objects in the globular 
cluster M22. They reported one clear microlensing event 
associated with an object of about 0.1 Mq, and six events 
unresolved in time that were very tentatively associated 
with planetary-mass objects in the cluster 

Observational evidence for mass segregation has been 
found in a number of globular clusters. In some cases, di- 
rect evidence comes from the measurement of a more cen- 
trally concentrated density profile for some heavier stellar 
population, such as blue stragglers or radio pulsars (Cote, 
Richer & Fahlman 1991; Layden et al. 1999; Rasio 2000). 
More often, the evidence comes from an apparent decrease 
in the slope of the (continuous) mass function toward the 
cluster core. This can either be estimated from the ob- 
served luminosity function at different radii (Howell et al. 
2001; Sosin 1997) or inferred from color gradients (How- 
ell, Guhathakurta, & Tan 2000). Note that, in contrast to 
younger star clusters where the present-day mass segrega- 
tion may still reflect initial conditions (Figer et al. 1999; 
but see Portegies Zwart et al. 2002), any mass segregation 
observed in the central regions of globular clusters must be 
a result of dynamical evolution, since the relaxation time 
in those regions is typically much shorter than globular 
cluster ages (by a factor ~ 10 — 10^). 

Our theoretical understanding of the dynamical evolu- 
tion of dense star clusters containing a wide mass spectrum 
is far from complete, even in the highly idealized limit of 
unevolving point masses. Direct A^-body simulations are 
very computationally intensive and have therefore been 
limited to systems with unrealistically low N ~ 10^ — 10* 
and rather narrow mass spectra (de la Fuente Marcos 
1996; Giersz & Heggie 1997; Takahashi & Portegies Zwart 
2000). Since realistic IMFs are thought to increase (per- 
haps steeply) toward smaller masses, a lower-mass cut-off 
(usually ~ 0.1 — 0.5 Mq) is introduced to avoid ending 
up with very few heavier objects and a large numerical 
noise in the simulation (since these heavier objects often 
dominate the overall dynamics). 

Instead, in this paper, we use Henon's Monte Carlo 
method to compute the dynamical evolution of clusters 
containing a more realistic number of stars (N = 10^). 
Our Monte Carlo code, as well as a number of test calcu- 
lations and comparisons with direct A/'-body integrations, 
have been described in detail in Joshi, Rasio, & Porte- 
gies Zwart (2000) and Joshi, Nave, & Rasio (2001). As 
a first step, we consider in this paper the simplest huo- 
component clusters, in which only two types of objects, 
of mass mi and m2 mi, are present. In addition, we 
assume that most objects in the cluster arc of mass mi, 
while the objects of mass m-i form a tracer population, 
i.e., the total component mass ratio M2/M1 <^ 1. In gen- 
eral, the gravothermal evolution of a two-component clus- 
ter can be either stable or unstable, depending on the two 
ratios m2/m,i and M2/M1 (Spitzer 1969). In the stable 
case, the two components remain thermally coupled and 
the cluster evolves dynamically toward energy equiparti- 
tion between the two species. In the unstable case, energy 



equipartition is impossible and the heavier species decou- 
ples thermally from the rest of the cluster, evolving sepa- 
rately toward core collapse and interacting with the lighter 
component through the mean gravitational field only. In a 
previous paper (Watters, Joshi, & Rasio 2000) we studied 
systematically the development of this instability (some- 
times called the Spitzer "mass-stratification" instability) 
in two-component systems, and we determined the loca- 
tion of the stability boundary in the parameter space of 
m2/m,i and M2/M1. Instead, in this paper, complemen- 
tary to our previous study, we concentrate on the stable 
systems with M2/M1 <^ 1 and we study systematically the 
effects of their dynamical evolution toward energy equipar- 
tition. In particular, we seek to characterize quantitatively 
the tendency for heavier objects to develop radial density 
profiles that are more centrally concentrated, as well as 
the tendency for lighter objects to concentrate away from 
the central regions, and to be preferentially evaporated in 
the Galactic tidal field. 

Our paper is organized as follows. In Sec. 2 we present 
an overview of our numerical approach, and we describe 
our initial conditions. In Sec. 3 we present our results for 
light tracers (with m2 < mi) and we discuss their impli- 
cations for microlensing by low-mass objects in globular 
clusters. In Sec. 4 we present our results for heavy tracers 
(with m2 > mi) and we compare them to simple analytical 
estimates. Our conclusions and summary are presented in 
Sec. 5. 

2. OVERVIEW OF NUMERICAL METHOD 

We study exclusively two-component clusters comprised 
of point masses, ignoring the effects of binaries and stel- 
lar evolution. To evolve the clusters, we use our recently 
developed 2-D Monte Carlo code (Joshi, Rasio, & Porte- 
gies Zwart 2000), which is based on Henon's Monte Carlo 
algorithm for solving the Fokker-Planck equation (Henon 
1971). Our code has shown close agreement with direct 
A-body and Fokker-Planck calculations of one-component 
clusters (Joshi, Rasio, & Portegies Zwart 2000), as well 
as calculations incorporating mass spectra and tidal mass 
loss (Joshi, Nave, & Rasio 2001). Furthermore, we have 
demonstrated reasonable agreement with the results of 
many different codes from Reggie's "Collaborative Exper- 
iment" (Heggie et al. (1998); see especially Fig. 2 of Rasio, 
Fregeau, & Joshi (2001)). 

For the initial conditions, we use the King model (a low- 
ered Maxwellian), given by the distribution function 



/ 



(e^/-^ - 1) 







for f > 
for f < 



(1) 



where £ = "i! — ^v^ is the (negative) energy of each star rel- 
ative to the potential at the tidal radius rt,"^ = —^+^{rt) 
is the relative gravitational potential, and k is a normaliza- 
tion constant. The quantity cr is a parameter, and should 
not be confused with the actual velocity dispersion y/Jv^- 
See, e.g., Binney & Tremaine (1987), or Spitzer (1987). 
For light tracers we enforce the tidal boundary induced by 
the Galaxy using a simple spherical Roche approximation 
(Joshi et al. 2001), while for heavy tracers (which tend 
to concentrate near the cluster center and hardly ever get 
ejected) we treat the clusters as isolated. 
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To effect a two-component distribution, we first create a 
single-component distribution of stars of mass mi accord- 
ing to equation (1), and then randomly replace a fraction 
of these background stars with tracers of mass m2, so that 
the initial density profile of the tracers is the same as that 
of the background. Although this creates clusters that are 
initially slightly out of virial equilibrium, we find that they 
settle into virial equilibrium within a few crossing times 
(a few timesteps). We have verified this independently us- 
ing direct iV-body integrations (sec below), which relax to 
virial equilibrium within a few crossing times. 

For light tracers (with mass ratio n = m2/m\ < 1), we 
consider a variety of King models for the initial conditions, 
with dimensionless central potentials Wa = 5'(0)/(j^ = 
{1, 3, 5, 7, 9}, a variety of mass ratios in the range 0.001 < 
H < 0.4, and a tracer population of either Nt = 1000 or 
Nt = 15000 out of 10^ total stars. Table 1 gives the initial 
conditions for each model. In all cases, tracers make up 
no more than 0.5% of the mass of each cluster, and conse- 
quently do not significantly affect the overall evolution of 
the cluster. 

For heavy tracers, we consider King models with Wq = 
{3, 7, 10}, a variety of mass ratios in the range 1.5 < ^ < 
10, and a tracer population in the range 100 < iVj < 1000. 
With a mass ratio greater than one, it is not always feasi- 
ble to use such a small number of tracers that the cluster 
is less than a few percent tracers by mass, since then sta- 
tistical noise is too great. Instead, we choose a reasonable 
value for the number of tracers, and consequently explore 
both systems which can reach energy equipartition, and 
systems which evolve away from thermal equilibrium. A 
simple test to determine whether thermal equilibrium is 
possible for a two-component cluster is the Spitzer stabil- 
ity condition (Spitzer 1969), 

3/2 

<0.16, (2) 



m2 
mi 



where Mi and M2 are the total masses of species 1 and 
2, respectively. If S is less than about 0.16, the system is 
"Spitzer-stable," and the two components will reach ther- 
mal equilibrium. If S exceeds this value, the system is said 
to be "Spitzer-unstable," and equipartition is impossible. 
Watters, Joshi, & Rasio (2000) performed a more refined 
analysis using our dynamical Monte Carlo code to evolve a 
wide range of two-component clusters, and arrived at the 
more accurate empirical stability criterion 

Table 2 gives these stability parameters, as well as the 
initial conditions, for each model considered. 

Some calculations with < 1 were repeated with di- 
rect A^-body integrations. For these direct A^-body cal- 
culations, we used the program kira within the Starlab 
software environment (see Portegies Zwart et al. (2001) 
and http : //manybody . org) and used the special- purpose 
GRAPE-4 and GRAPE-6 hardware (a single GRAPE-6 
board) (Makino et al. 1997) to accelerate the computation 
of gravitational forces between stars. Time integration of 
stellar orbits was accomplished using a fourth-order Her- 
mite scheme (Makino & Aarseth 1992). kira also incor- 
porates block timesteps (McMillan 1986a,b; Makino 1991) 
and a special treatment of close two-body and multiple 



encounters of arbitrary complexity. Given the high cost 
of the direct A^-body integrations we limited our selec- 
tion of initial conditions to a few models with /x < 1 and 
performed the calculations with N = 6144 stars; these cal- 
culations take of order 10 days each to complete on either 
the GRAPE-4 or a single GRAPE-6 board. (With our 
Monte-Carlo code, each calculation takes ^ 12 hours.) 

Following the convention of most previous studies, we 
define dynamical units so that G = Mq = — 4i?o = 1, 
where Mq and E{] are the; initial total mass and total en- 
ergy of the cluster (Henon 1971). The units of length, L, 
and time, T, are then 

i = GMo^(-4£;o)-S and T = GmI''^{-AEo)-^I'^ . 

Application of the virial theorem shows that L is the virial 
radius of the cluster and T is of the order of the initial dy- 
namical (crossing) time. In the code, we use as the unit of 
time the initial relaxation time, tr, which is given by 

i. = T%T, (5) 
m Ao 

where A^o is the initial total number of stars in the cluster. 
For a two-component cluster, it is possible to define sim- 
ilar relaxation times for each component separately. The 
relaxation time for component 1, tri, is given by 

,6) 



trl 



\nN^ 



where A'"i is the initial number of stars of component 1, and 
Ti is the crossing time for component 1 , given by equation 
(4) with Afo and Eq replaced by the initial total mass and 
energy of component 1, respectively. The relaxation time 
for component 2,tr2, is given analogously, although the re- 
sulting expression is only applicable when the two species 
have decoupled (e.g., a collection of black holes undergoing 
the Spitzer mass-stratification instability). 

To facilitate comparison with previous studies, we re- 
port times in this paper in units of the initial half-mass 
relaxation time, trh, given by the standard expression (see, 
e.g., Spitzer 1987) 



trh = 0.138- 



7\rl/2 3/2 



(7) 



mV2GV2lnAro ' 
where rh is the initial half- mass radius, and m = Mq/Nq 
is the average stellar mass. 

Each calculation is terminated when the 0.35% La- 
grangian radius of the heavier component falls below 0.001 
(in our length units). The Lagrangian radii are then in- 
spected graphically and the core collapse time, tec, de- 
termined by noting the time at which the innermost La- 
grangian radii of the heavier component begin to dip ap- 
preciably. Core collapse times determined in this way are 
not sensitive to the precise criterion used but can have a 
large statistical uncertainty, particularly when the number 
of stars within the innermost Lagrangian radii is small. 

3. LIGHT TRACERS 

Intuition has it that globular clusters cannot contain 
extremely light stars, since equipartition of energy would 
imply a mean velocity for the light stars that far exceeds 
the escape velocity for the cluster (Taillet, Longaretti, & 
Salati 1995). Although proportionately more light tracers 
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Fig. 1. — Comparison of Monte Carlo and Af-body data. (Solid line is Monte Carlo, dotted line is Af-body.) Shown are the 5, 10, 25, 
50, 75, and 90% Lagrangian radii of background stars (left plot) and tracers (right plot) for a Wo = 5 King initial model with a mass ratio 
/i = 0.001, showing reasonable agreement between the two methods. 



than background stars are lost during the evolution, we 
find that a significant fraction arc retained in the halo. In 
some cases their number density there actually increases 
during the evolution. 

3.1. Numerical Results 



To ensure the validity of our Monte Carlo results, we 
first sought to compare integrations performed with the 
Monte Carlo code with the results of direct A'^-body sim- 
ulations. As a representative case, Figure 1 shows 5, 10, 
25, 50, 75, and 90% Lagrange radii (radii containing a 
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Fig. 2a. — Lagrangian radii of background stars (solid line) and 
tracers (dotted line) from Figure 1, showing only the Monte Carlo 
results. Note the stark contrast between the two, illustrating the 
extreme mass segregation present even early in the simulation, at 
t ~ 5trh- By core collapse, over 90% of the tracers occupy the 
outer region of the cluster in which 50% of the background stars 
reside. 



Fig. 2b. — Same as Figure 2a, but for a very centrally con- 
centrated Wo = 9 King initial model. Even though this cluster 
reaches core collapse very quickly (because it is initially almost 
core-collapsed), by that time over 90% of the tracers occupy the 
outer region of the cluster in which 75% of the background stars 
reside. 
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Fig. 2c. — Same as Figure 2a, but for a mass ratio // = 0.1. Al- 
most the same degree of mass segregation present for the extreme 
mass ratio of fi = 0.001 is witnessed for this more reasonable value. 
Again, by core collapse, about 90% of the tracers occupy the outer 
region of the cluster in which 50% of the background stars reside. 



Fig. 3a. — Number histogram of stars versus radius for a Wo = 5 
King initial model with a mass ratio fj, = 0.001, showing back- 
ground stars in the top plot and tracers in the bottom. The dotted 
line represents a time very near the start of the simulation, while 
the solid line represents a time midway through the simulation, at 
t = 4.9 (rh- While this isn't a number density plot, it gives a feel 
for the radial profile of the two species. 



constant mass fraction) for background stars (left) and 
tracers (right) for a Wo = 5 King initial model with mass 
ratio fi = 0.001. The agreement remains quite good at all 
times. The two methods agree remarkably in the most rel- 
evant result: that the light tracer stars "diffuse" into the 
outer halo, while the core, comprised mainly of heavier 
"background" stars, contracts. This is consistent with the 
mass segregation found with previous A^-body and Fokker- 
Planck studies (Spitzer 1987). 

Figure 2a displays the segregation more clearly, showing 
only the Monte Carlo results from Figure 1. Here back- 
ground stars are represented by a solid line, light trac- 
ers by a dotted line. As one expects from simple energy 
equipartition arguments, the tracers gain energy from in- 
teractions and are pushed out of the core, as is seen in the 
mass segregation evident as early as t ~ 5trh- However, 
25% of the tracer population is left to linger in the outer 
regions at the time of core collapse. (By contrast, 65% of 
the background stars remain in the cluster.) Instead of 
receiving large monolithic kicks from the heavy stars and 
being immediately cjec;ted from the cluster, the light trac- 
ers slowly gain energy and are gradually pushed out into 
the halo, where there are so few heavy stars that the light 
tracers are unable to exchange energy and are left there 
to linger. Figure 2b is the same as Figiire 2a, but for a 
very centrally concentrated Wq — 9 King initial model. 
Because this cluster is initially so close to core collapse, it 
does not have much time to evolve; yet by collapse over 
90% of the tracers occupy the outer region of the cluster in 
which 75% of the background stars reside. This extreme 
segregation is witnessed even for a more moderate mass 
ratio of /U = 0.1. Figure 2c shows the same as Figure 2a, 



but for a Wq = 5, /x = 0.1 model. Here, by the time of 

core collapse, about 90% of the tracers occupy the outer 
region in which 50% of the tracers reside. 

For a more detailed look at the radial profile of the light 
tracers relative to the background, Figure 3a shows the 
number of stars in each radial bin. The number in each 
bin is divided by the initial total number of stars in the 
cluster, iVo, so that the vertical axis represents the fraction 
of stars in each bin. This histogram, which is proportional 
to is not a number density plot but still gives a feel 
for the number density, and was chosen because it most 
clearly displays phenomena at large radius. Shown is the 
same Wq = 5, /z = 0.001 model as in Figure 1, with back- 
ground stars in the top plot and tracers in the bottom. 
The dotted line represents a time very near the start of 
the simulation, while the solid line represents a time mid- 
way through the simulation, at t = 4.9 t^h- There is a clear 
increase in the density of background stars in the core dur- 
ing the evolution, while there is also a clear decrease for 
tracers. More importantly, there is a significant enhance- 
ment of tracers in the halo region. Similarly, Figure 3b 
shows the same results for a Wq = 7, /i = 0.001 model. 
In this case, the number density of tracers in the halo is 
very visibly greater at core collapse than at the start of 
the simulation, and by about 50% in some regions. Even 
for the more modest mass ratio of /i = 0.1, a similar profile 
is obtained, as shown in Figure 3c. Although a bit noisier 
than the preceding, since only 1000 tracers are used here, 
the same significant enhancement of tracers in the halo is 
evident. 

Table 1 gives relevant initial conditions for all the light 
tracer systems considered here, as well as core collapse 
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Fig. 3b. — Same as Figure 3a, but for a Wo = 7 King ini- 
tial model. Here the solid line represents a time near collapse, 
t = 7.3 trh- In this case, the number density of tracers in the halo 
is very visibly greater at core collapse than at the start of the 
simulation, and by about 50% in some regions. 



Fig. 3c. — Same as Figure 3a, but for a Wo = 7 King initial 
model with a mass ratio n = 0.1. Again, the solid line represents 
a time near core collapse, t = 7.2 4^)^. Although the curve is a bit 
noisy, there is clearly an enhancement in the number density of 
tracers in the halo, even at this modest mass ratio. 



times, the fraction of each species left in the cluster at 

core collapse, and the ratio of half-mass radii for the two 
species at core collapse. Wc have also included a few as- 
trophysically relevant timcscalcs: Ti.5(rc), the time for the 
number of light stars within the initial core radius to de- 
crease by a factor of 1.5; Tio(rc), the time for the number 
of light stars within the initial core radius to decrease by a 
factor of 10; Ti.5(r/i), the time for the number of light stars 
within the initial half-mass radius to decrease by a factor 
of 1.5; and rio(r/i), the time for the number of light stars 
within the initial half-mass radius to decrease by a fac- 
tor of 10. Figure 4 shows a representative determination 
of these timcscalcs. Smooth curves were fit to the data 
using gnuplot's 'acsphnes' (approximation cubic splines) 
routine. The timescales show a clear trend, increasing as 
the mass ratio approaches unity, consistent with energy 
equipartition arguments. A closer look reveals a critical 
mass ratio, around 0.3, below which the timescales no, 
and especially T1.5, are roughly constant: the light stars 
are so light that they are immediately ejected from the in- 
ner regions of the cluster on a timescale that is independent 
of mass ratio. A simple derivation of this critical mass ra- 
tio, /Ucrit, proceeds as follows. The escape speed from the 
center of the cluster is given by v'^ = 2^'(0) = 2Wo<t'^- 
In the core of the cluster, densities are very high and it 
is expected that the two stellar species will attain ther- 
mal equilibrium: (w^)core = M(t'|)core- Since the cluster 
is dominated by the heavy species (species 1), its velocity 
dispersion can be obtained directly from the distribution 
function (equation (1)): 



Jo 



V2*(0) 



(8) 



This equation can be integrated numerically to give 
('fi)coro in units of cr^: (wx)coro = acr^. Setting the ve- 
locity dispersion of the lighter species equal to the escape 
speed at the center yields a simple equation for the critical 
mass ratio, 

a , , 

" 2Wo ' ^ ' 

below which the light stars are almost immediately ejected 
from the core (within a core relaxation time) , but not nec- 
essarily from the cluster. Table 3 gives a and /icrit for 
many values of Wq. Comparison with Table 1 shows that, 
indeed, below this mass ratio the light stars arc ejected 
from the central regions of the cluster on a timescale that 
is independent of mass ratio; above this mass ratio, the 
segregation timescale increases dramatically. 

3.2. Implications for microlensing 

It is possible that globular clusters could contain co- 
pious amounts of non-luminous matter often referred to 
as (baryonic) dark matter in this context (Heggie, Griest, 
& Hut 1993; Taillet, Longaretti, & Salati 1995; Heggie & 
Hut 1996). However, it is probable that the amount of 
dark matter is poorly constrained by observations of lumi- 
nous stars. Consider a cluster with a continuous Salpeter 
IMF, 

4S oc . (10) 

Since we want to study very light objects, let the IMF ex- 
tend from 0.001 Mq to 10 Mq (even though such a mass 
function was never meant to apply below ~ O.IMq!). To 
place this problem in the context of two-component clus- 
ters, break up the IMF into two bins: a tracer bin ex- 
tending from 0.001 Mq to 0.1 Mq corresponding to non- 
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Table 1 



Main results for light tracers 



Wo 




Nt 


Nt 


'cc 
'rh 


"^l.sC^c) 








JVtCcc) 


NbCco) 


'■htCco) 


'rh 


'rh 


'rh 


'rh 






'■hhCco) 


1 


0.001 


1000 


99000 


9.7 


1.54 


9.25 


1.83 


8.44 


0.039 


0.32 


5.6 




0.001 


15000 


85000 


8.5 


1.67 


8.40 


1.41 


7.63 


0.044 


0.32 


5.5 




0.004 


1000 


99000 


9.7 


1.92 


8.90 


2.05 


8.48 


0.040 


0.33 


4.4 




0.01 


1000 


99000 


9.4 


1.70 


*** 


1.71 


8.64 


0.051 


0.33 


4.1 




0.04 


1000 


99000 


9.6 


1.65 


8.99 


1.73 


8.64 


0.048 


0.32 


4.6 




0.1 


1000 


99000 


9.6 


2.45 




2.41 




0.055 


0.33 


3.6 




0.2 


1000 


99000 


9.6 


2.64 


*** 


2.22 


9.35 


0.055 


0.32 


3.2 




0.4 


1000 


99000 


9.6 


6.80 


*** 


4.49 


*** 


0.093 


0.33 


2.7 


3 


0.001 


1000 


99000 


11.7 


0.90 


11.63 


1.44 


10.45 


0.076 


0.46 


2.9 




0.001 


15000 


85000 


10.1 


0.87 


*** 


0.97 


8.70 


0.086 


0.46 


3.2 




0.004 


1000 


99000 


11.6 


1.05 


11.13 


1.41 


10.57 


0.079 


0.46 


2.9 




0.01 


1000 


99000 


11.7 


1.08 


10.63 


1.65 


10.13 


0.068 


0.46 


3.2 




0.01 


15000 


85000 


10.0 


1.13 




1.03 


8.90 


0.088 


0.46 


3.0 




0.04 


1000 


99000 


11.6 


1.46 




1.55 


10.93 


0.10 


0.46 


3.1 




0.042 


15000 


85000 


10.2 


1.09 


*** 


1.13 


9.37 


0.094 


0.46 


3.0 




0.1 


1000 


99000 


11.7 


1.05 


*** 


1.74 


*** 


0.11 


0.46 


2.8 




0.2 


1000 


99000 


11.7 


1.99 




2.35 


11.65 


0.11 


0.46 


2.5 




0.4 


1000 


99000 


11.7 


10.61 


*** 


5.61 




0.16 


0.47 


2.2 


5 


0.001 


1000 


99000 


13.3 


0.43 


11.4 


1.05 


11.5 


0.26 


0.64 


2.8 




0.001 


15000 


85000 


11.1 


0.38 


10.9 


0.79 


9.54 


0.27 


0.65 


2.8 




0.1 


1000 


99000 


12.9 


1.43 


*** 


1.52 


12.78 


0.30 


0.65 


2.5 


7 


0.001 


1000 


99000 


8.7 


0.57 


8.31 


0.87 


7.22 


0.68 


0.87 


2.3 




0.001 


15000 


85000 


7.6 


0.21 


2.30 


0.70 


7.25 


0.68 


0.87 


2.2 




0.004 


1000 


99000 


9.0 


0.67 


7.24 


1.07 


8.20 


0.67 


0.86 


2.3 




0.01 


1000 


99000 


8.9 


0.50 


6.90 


1.00 


8.77 


0.68 


0.87 


2.1 




0.04 


1000 


99000 


9.0 


0.54 


6.20 


1.08 


*** 


0.68 


0.87 


2.2 




0.1 


1000 


99000 


8.9 


0.54 


8.44 


1.17 




0.70 


0.87 


2.0 




0.2 


1000 


99000 


8.9 


0.70 


*** 


1.28 


+ * + 


0.73 


0.87 


2.1 




0.4 


1000 


99000 


8.8 


1.76 


*** 


2.51 


*** 


0.77 


0.87 


1.8 


9 


0.001 


15000 


85000 


1.1 


0.019 


0.165 


0.62 


*** 


0.97 


0.98 


1.4 




0.1 


1000 


99000 


1.5 


0.095 


0.875 


1.02 


*** 


0.97 


0.98 


1.4 



Note. — Here Nt is the number of trace stars; A'^;, is the number of background stars; tec is the core collapse time; ri.6(r) 
is the time it takes for the number of tracers within radius r to decrease by a factor 1.5; rio(r), by 10; and rht and rhb are 
the half-mass radii of the tracers and background stars, respectively. represents data that were not attained before core 

collapse. 



Table 2 



Main results for heavy tracers 



Wo 


M 


Nt 


Nt 


S 


A 


'cc 
'rh 


^ol.'-h<'-<'-0 


'rh 


'rh 


-1.3(-h) 
'rh 


-htCcc) 
-hhCcc) 


3 


1.5 


1000 


99000 


0.028 


0.040 


16.50 


0.74 


0.52 


11.61 


1.30 


0.070 




3 


1000 


99000 


0.16 


0.42 


7.60 


0.27 


0.33 


3.85 


0.54 


0.099 




6 


100 


99900 


0.088 


0.44 


8.30 


0.19 


0.25 


1.34 


0.39 


0.095 




10 


100 


99900 


0.32 


2.5 


3.20 


0.13 


0.42 


1.61 


0.44 


0.080 




10 


1000 


99000 


3.2 


25.4 


1.50 


0.16 


0.37 


1.36 


0.37 


0.18 


7 


1.5 


1000 


99000 


0.028 


0.040 


8.80 


0.54 


0.31 


5.16 


3.91 


0.23 




3 


1000 


99000 


0.16 


0.42 


3.10 


0.22 


0.053 


0.81 


0.84 


0.098 




6 


1000 


99000 


0.89 


4.5 


0.50 


0.13 


0.054 


0.37 


0.39 


0.54 




10 


1000 


99000 


3.2 


25.4 


0.31 


0.09 


0.032 


0.25 


0.27 


0.60 


10 


1.5 


200 


99800 


0.0055 


0.0080 


0.38 


0.31 


0.028 


0.37 




0.85 




1.5 


1000 


99000 


0.028 


0.040 


0.38 


1.5 


0.028 


0.32 


*** 


0.98 




3 


200 


99800 


0.031 


0.084 


0.19 


0.17 


0.0021 


0.043 


*** 


0.86 




3 


1000 


99000 


0.16 


0.42 


0.081 


0.6 


0.0021 


0.042 




0.98 




6 


200 


99800 


0.18 


0.89 


0.022 


0.09 


0.0014 


0.012 




0.92 




6 


1000 


99000 


0.89 


4.5 


0.013 


0.4 


0.0016 






1.0 




10 


100 


99900 


0.32 


2.5 


0.021 


0.04 


0.00050 


0.0055 


*** 


0.91 



Note. — Variables are the same as Table 1, except that now ri.sfr) is the time it takes for the number of tracers within 
radius r to increase by 1.5. Also included are the mass segregation timescale tq (in units of tr), as well as the stability 
parameters S and A. Note that for each simulation the number of stars lost from the cluster was 1% or less. 
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Fig. 4. — A representative determination of the mass segregation 
timescales Ti.6(rh) and Tio(rh) from Table 1 for the case /n = 0.004, 
Wo = 3. The top horizontal line is the number of trax;ers initially 
within the half-mass radius, No{rh); the middle horizontal line 
Noi'rh)/^-^, and the bottom horizontal line Afo(r;i)/10. 



Table 3 

Velocity dispersion and critical mass ratio 



VVQ 


<.?>ccre 


Merit 


1 


0.799 


0.399 


2 


1.466 


0.367 


3 


1.990 


0.332 


4 


2.373 


0.297 


B 


2.632 


0.263 


6 


2.796 


0.233 


7 


2.892 


0.207 


8 


2.945 


0.184 


9 


2.973 


0.165 


10 


2.987 


0.149 


11 


2.994 


0.136 


12 


2.997 


0.125 



Note. — Velocity disper- 
sion in the core in units of 
the velocity dispersion pa- 
rameter (T, and critical mass 
ratio for immediate ejection 
from the core for many val- 
ues of the dimensionless cen- 
tral potential Wq; see the 
discussion preceding equa- 
tion (9). 



luminous objects, and a background bin extending from 
0.1 Mq to 10 Mq corresponding to main sequence stars. 
(For simplicity, wc will ignore stellar evolution and sup- 
pose that the heavy stars keep their mass for the lifetime 
of the cluster.) Integrating equation (10) over these two 
bins gives a mean tracer mass of 0.0031 Mq, a mean back- 
ground mass of 0.31 Mq, and a ratio of the number of 
tracers to background of Nt/Ni, = 500. Using the back- 
ground stars as species 2 in equations (2) and (3), we find 
S = 200 and A = 1.3 x 10^: this simple model is well 
outside the Spitzer stability regime. The two species will 
decouple completely and evolve independently; observa- 
tions of the luminous background stars will not betray the 
presence of the very light tracers. 

Approximating the continuous mass spectrum by two 
bins is most likely unjustified, however, since in a continu- 
ous mass spectrum stars of mass M will certainly exchange 
energy with stars of mass M + dM . Thus light stars will 
exchange energy with slightly heavier stars, which will ex- 
change energy with slightly heavier stars, etc., providing 
an indirect means for light objects such as Jupiters to ex- 
change energy with heavy objects such as main-sequence 
stars. Thus thermal equilibrium among all species may 
not be forbidden in a realistic cluster, but the time scale 
over which it occurs is unclear. 

Taillet, Longarctti, & Salati (1995) and Taillet, Salati, & 
Longaretti (1996) considered multi-component King mod- 
els and assumed that thermal equilibrium did occur among 
species. Using analytical methods, they found that the 
cluster surface brightness was nearly unaffected by the 
amount of non- luminous stars present, for a cluster with 
dark matter content ranging from 0% to 90% by mass. 
Therefore, observations of the luminous cluster stars would 
give no evidence of the dark matter in the cluster. (How- 
ever, in the most extreme cases, the projected velocity 
dispersion profiles could be used to detect dark matter by 
observing the luminous component.) 

Regardless of whether thermal equilibrium occurs in 
clusters, it is probable that dark matter would not be eas- 



ily detectable by observing the luminous members of clus- 
ters. Dark matter could be directly detected, however, by 

its gravitational lensing effect on distant, luminous sources 
(Paczyhski 1986, 1991, 1994). By convention, a microlens- 
ing event is said to occur when a background source lies, 
in projection, within the Einstein radius of a lens and is 
consequently amplified by a factor greater than 1.34 (Refs- 
dal 1964; Victri & Ostrikcr 1983). The optical depth, t, is 
the probability that such an event will occur, and for low 
probability, is given by 

T = ^ JdVn{Dd)nel , (11) 



where 5u) is the solid angle of sky considered, dV = 
SujD^dDd is the clement of volume, n is the number den- 
sity of lenses, Dd is the distance to the lens (defiector), and 
9e is the angular Einstein radius of the lens. The Einstein 
radius is given by the usual expression 
2 4GM - Dd 
D,D, ' 

where Dg is the distance to the source. Here wc will be 
considering lensing of sources in the Galactic bulge, for 
example, by globular clusters along the line of sight to the 
bulge; therefore, typical values for Dd and Dg are ~ 3 kpc 
and ~ 8 kpc, respectively. Since tidal radii of globulars 
arc typically of order tens of parsccs, wc can approximate 
Dd and Dg as constant over the range of integration to 
arrive at the simpler expression 



47rG 



T = 



S(r)Dd 1 - 



Dd 
Ds 



(13) 



where S(r) is the surface density along the line of sight, 
and r is the distance perpendicular to the line of sight. 

Figure 5a shows the optical depth due to heavy stars 
(solid line) and light stars (dotted line) for a Wq = 3 ini- 
tial King model with 15000 tracers of mass ratio fi = 0.01. 
We have assigned the background stars a mass of 0.31 
Mq, and multiplied the tracer number density by a factor 
500 X (85000/15000), consistent with the preceding discus- 
sion of approximating the Salpeter IMF by two mass bins. 
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Fig. 5a. — Optical depth to microlenslng due to heavy stars (solid 
line) and light stars (dotted line) for a Wo = 3 initial King model 
with 15000 tracers of mass ratio /x = 0.01. The number density 
of light stars has been amplified to emulate a continuous Salpeter 
IMF. This snapshot was taken near core collapse, at t = 9.5 trh- 
Generic values of D^i = 3kpc and Ds = 8kpc were used in the 
calculation. 
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Fk;. 5d. Optical depth to microlenslng at present due to lumi- 
nous stars (solid line) and dark matter (dotted line) for the globular 
cluster M22. 



(This heuristic model therefore defines a thermally equilib- 
rized two-component idealization of the continuous mass 
function.) We have also multiplied both optical depths by 
a factor of 5 so that the figure is normalized to a cluster 
that initially had 5 x 10^ stars and a mass of 1.3 x 10^ Mq. 
Generic values of Da = 3 kpc and Dg = 8 kpc were used in 
the calculation, and this snapshot of the system was taken 
at t = 9.5 trh- Note not only the mass segregation, but 
also that at large radii the optical depth duo to the light 
stars is about an order of magnitude greater than that due 
to the heavy stars. Significant numbers of non- luminous 
stars may linger in the outer regions of globular clusters. 

3.3. Dark matter in M22 

Sahu et al. (2001) have carried out ffST observations of 
the central regions of the globular cluster M22 (NGC 6656) 
and have detected at least one microlenslng event. All the 
parameters needed to deduce the mass of the lens from the 
timescale of the lensing event are fairly well known so that 
a prediction of its mass can be made. Assuming distances 
of 8.2 kpc to the source (for a typical star in the Galactic 
bulge) and 2.6 kpc to M22, as well as a proper motion for 
M22 of 10.9 mas yr"^ (134.4 km s^^), they determine the 
mass of the lens to be O.ISIq qjM©. This result is exciting 
because it is the first such detection where the mass of the 
lens can be estimated. They also report a tentative detec- 
tion of 6 sub-Jupiter mass objects in the core (an upper 
limit of ~ 0.25 times the mass of Jupiter for each) that 
has since proven spurious (Sahu, Anderson, & King 2001). 
Although the events are no longer associated with planets, 
we still present our theoretical analysis of the amount of 
dark matter in M22 as a useful exercise, illustrating how 



future detections could be used to estimate the amount of 
dark matter in clusters. 

It is a simple matter to calculate the optical depth due 
to dark matter given these detections. Since the lenses are 
moving with respect to the source, the actual probability 
for microlenslng is modified from equation (11) to 



Prob = = — 

No Slu 



dVn{Dd)29EipAt, (14) 



where is the number of microlenslng events observed, 
No is the number of observations, tp is the angular velocity 
of the lens, and At is the time duration of each observa- 
tion. In terms of the optical depth, this is simply 

- - (15) 



No 



2 At 

TT to 



(16) 



where to = OE/<f is the characteristic timescale of a mi- 
crolenslng event. The optical depth is then 

^ nNeto_ 
2 No At' 

Sahu et al. (2001) monitored approximately 83,000 stars 
for a period of ~ 105 days, and found 6 candidate dark 
matter (below 0.1 Mq) microlenslng events lasting at most 
0.8 days each. Equation (16) then gives an optical depth 
due to dark matter of 8.7 x 10^^. Adopting the same 
simple prescription used before of approximating a con- 
tinuous mass function by two mass bins, we can evolve a 
two-component cluster from a reasonable initial state until 
its profile resembles that of M22. We can then adjust the 
number of tracers so that the optical depth in the core due 
to dark matter matches observations, and consequently es- 
timate what fraction of the cluster mass would have been 



10 



FREGEAU ET AL. 



in dark matter had the 6 tentative detections proven cor- 
rect. 

For our initial conditions we take a Wq = 3 King model 
with 85000 background stars of mass 0.38 Mq and 15000 
tracers of mass 0.016 Mq to emulate a Kroupa IMF — 
which we have extended well below the 0.1 Mq lower limit 
for which it was intended — from 0.001 Mq to 10 Mq split 
into two bins at 0.1 Mq (Kroupa, Tout, & Gilmore 1993). 
The Kroupa IMF is more conservative than a Salpeter in 
that it contains far fewer sub-stellar-mass objects. (The 
exact form of the IMF chosen is largely irrelevant, since 
in the end we scale the number of stars in each mass bin 
to match observations; however, it may affect the percent- 
age of light stars lost from the cluster during evolution.) 
We then evolve the cluster until its profile matches that 
of M22, approximately a Wq = 6 model (Harris 1996). 
We adjust the number of background stars so that the 
mass due to luminous matter is 10'''"'''''AfQ (Mandushev, 
Spassova, & Staneva 1991), and adjust the number of 
tracer stars so that the optical depth due to dark mat- 
ter in the core is approximately 8.7 x 10~^. 

Figure 5b shows the optical depth due to luminous stars 
(solid line) and dark matter (dotted line) using the pre- 
scription just given. Adopting the total cluster mass and 
dark matter optical depth quoted above, we find that the 
current mass of M22 would have been ~ 20% dark matter 
by mass. If either the total cluster mass or optical depth 
are off by a factor of two, this number could be as small 
as 12% or as high as 35%. If both are off by a factor of 
two, then this number could be as small as 6% or as high 
as 52%. This range of values is consistent with the value 
quoted by Sahu et al. (2001) of ^ 10%, but implies a higher 
median. Using the median values, this analysis implies 
that M22 would have had an initial mass of 8.3 x 10^ Mq, 
35% of which was dark matter (this figure could be any- 
where between 12% and 69%), and that it would have lost 
2 X 10^ Mq of dark matter to the Galaxy. These ranges of 
values should be treated with caution, however, since our 
Monte-Carlo code ignores the effects of binaries, mass loss 
due to stellar evolution, and the effects of a non-spherical 
tidal field, all of which may affect the fraction of light trac- 
ers retained in a cluster. Furthermore, it is clear that in 
order to tighten constraints on the amount of dark matter, 
one must also refine estimates of the distances to the lens 
and the source. 

3.4. Upper Limit on Dark Matter in M22 

All six "planetary" events are no longer interpreted as 
being due to microlensing. More careful examination re- 
veals that the apparent brightening observed in both of a 
pair of cosmic-ray-split images was actually due to point- 
like cosmic rays hitting close to the same star in both im- 
ages (Sahu, Anderson, & King 2001). This null detection, 
while unfortunate, can be used to estimate an upper limit 
on the mass fraction of M22 in very low-mass objects. 

Adopting our simple prescription of approximating the 
continuous mass function by two bins, one containing lu- 
minous matter (objects heavier than 0.1 Mq), the other 
containing dark matter (objects lighter than 0.1 Mq), we 
assume that the dark matter in M22 is represented by a 
collection of objects with the average mass in that bin, 
0.016 Mq. We then expect the number of dark matter 



lensing events in an observation to obey a Poisson distri- 
bution. Adopting a mean of 1 dark matter lensing event, 
using ~ 6 days for a 0.016 Mq lens, and performing es- 
sentially the same analysis as in Sec. 3.3, a null detection 
yields an upper limit on the current mass fraction of M22 
in very low-mass objects of ~ 25% at the 63% confidence 
level. Adopting a mean of 3 events yields an upper limit 
of ~ 50% at the 95% confidence level. 

Although our simple calculation is physically motivated, 
it is necessarily crude, sinc;e a (;arc!ful analysis should take 
into account the mass distribution of the dark matter. In 
addition, there is another effect to consider. The typical 
separation between exposures in the Sahu et al. (2001) ob- 
servation is about 3 days. Thus a microlensing event with 
characteristic width smaller than 3 days, corresponding to 
a lens mass of 0.004 Mq, has a non-zero probability of not 
being detected, biasing detections toward higher masses. 
Since the smallest mass we consider is 0.001 Mq, presum- 
ably this effect is small. 

4. HEAVY TRACERS 

4.1. Analytic Results 

Spitzer (1969) pioneered the theoretical analysis of mass 
segregation in clusters. Considering a two-component clus- 
ter in which one component is more massive than the 
other, yet sparse enough to be considered a tracer pop- 
ulation, he found that the timescale on which the heavier 
species contracts varies approximately as 1//Lt. For com- 
pleteness, we include a brief summary of the steps leading 
to this result. 

For the generic case of energy transfer when E2, the 
mean kinetic energy of the massive stars, is greater than 
the corresponding energy Ei for the lighter stars, the rate 
of energy loss of the massive stars is given by 

dU2 _ {E^ - E2) ^ ^^^^ 



dt 



t. 



eq 



where U2 is the mean total energy of the heavier species. 
The equipartition time is approximately (Spitzer 1940, 
1962) 



3/2 



8 (67r)^/ Vol G^ma IniVi 



(18) 



where poi is the central density of the lighter species, A''i 
is the number of light stars in the cluster, and brackets 
denote an average over stars. When one expresses this 
timescale in units of t^i , the overall relaxation timescale of 
species 1, one obtains 



t, 



_eq 



37rV2 



16 m2 V 

It follows that if the rms velocities of the two species are 
equal, as they are at the start of our simulations (they re- 
main roughly equal for the first half of the simulation; see 
Figure 7) , teq/iri initially varies approximately as l//x. 

To see how this timescales manifests itself, assume, for 
the sake of simplicity, that the kinetic energy of the heav- 
ier species, E2, dominates that of the lighter species, Ei. 
Equation (17) then becomes 

(20) 



dt 



t, 



eq 
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Fig. 6. — Number of tracers in the halo {r^ < r < rt) as the 
cluster contracts, showing a representative determination of the 
mass-segregation timescale tq for the case Wq = 3, h = 1.5, and 
Nt = 1000. 
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Fig. 7. — Ratio of rms velocities of the heavier component to the 
lighter component for the system shown in Figure 6. Solid line is 
this ratio for the cluster as a whole. Note that this ratio remains 
near unity for about the first half of the simulation. For reference, 
the value of this ratio in the core is shown as a dotted line. This 
curve approaches 2/3, consistent with equipajtition. 



Application of the usual virial theorem yields U2 = —E2 = 
\W2, where W2 is the potential energy of species 2, and 
so we have 



W2 oc e*/*«i . 



(21) 



Under the assumption that the potential of species 2 dom- 
inates, W2 is proportional to GM2/r2h, where r2h repre- 
sents any characteristic radius of species 2, but for con- 
creteness has been chosen to be the half-mass radius. 



0.1 



' 1 ' 


1 


- „Wo=3 




^ \*o = '?'\\ 

- X; 




w„=io\^^ 








! 0.8//J. 




0.5/;ti 

, 1 





1 


10 




t/t. 



Fig. 8. — Mass-segregation timescale as a function of mass ratio, 
for a few King initial models. In all cases to oc as predicted 
by simple theoretical arguments. 



Fig. 9. — Ratio of half-mass radii of the tracers to the background 
stars for the system shown in Figure 6, displaying persistent segre- 
gation. Also plotted, for reference, is the overall half-mass radius 
of the cluster. 
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Thus, 

r2h oc e-*/*«"> , (22) 

which simply states that the timescalc for contraction of 
the heavier species is given by tgq, the cquipartition time. 
Equation (22) still holds true even when the potential of 
species 2 is not the dominant potential. 

4.2. Numerical Results 

For studying mass segregation of heavy tracers, wc use 
as the unit of time the overall cluster relaxation time, tr, 
given by equation (5). This time differs from tri (given 
by equation 6), the relaxation time of species 1 and the 
unit of time in equation (19), by a factor close to unity 
that varies slowly among sets of initial conditions. The 
quantity tr also differs from the half-mass relaxation time, 
trh (given by equation 7), the unit of time used in the pre- 
ceding sections for studying light tracers. A good rule of 
thumb is that trh ~ O.ltr (Joshi, Rasio, & Portegies Zwart 
2000). 

To determine the mass segregation timescale — which 
should approximately equal icq from equation (22), but 
which wc label tq wc look at the number of tracers in 
the halo, defined to be the region between the half-mass 
radius and the tidal radius, as a function of time and fit to 
it a decaying exponential. Figure 6 shows a representative 
determination of this timescalc for a Wq = 3, fi = 1.5, 
Nt = 1000 model. Note that the exponential fit is only 
meaningful for the first few half- mass relaxation times, 
during which the ratio of velocity dispersions (^^/(I'f) 
doesn't change appreciably. Figure 7 shows this ratio hov- 
ering at unity for nearly the first half of the simulation. 
(For reference, {v2)/{vi) averaged over the core is shown 
as a dotted line. This curve approaches 2/3, consistent 
with thermal equilibrium.) Figure 8 shows tq as a func- 
tion of /i for King initial models with Wq = 3, 7, and 10. 
In all cases, the mass segregation timescale shows a clear 
dependence, as expected from equation (19). Figure 
9 gives another look at the segregation, showing the ratio 
of half-mass radii of the tracers to the background stars. 

Table 2 gives relevant initial conditions for all the heavy 
tracer systems considered here, as well as the thermal equi- 
librium stability parameters S and A (see eqs. (2) and (3)), 
core collapse times, and the ratio of half-mass radii for the 
two species at core collapse. Similarly to Table 1 we have 
included the following astrophysically relevant timcscalcs: 
Toifh < r < rt), the e-folding timescale shown in Figure 8; 
Ti.5{fc), the time for the number of heavy stars within the 
initial core radius to increase by a factor of 1.5; Tio{rc), the 
time for the number of heavy stars within the initial core 
radius to increase by a factor of 10; and Ti.5(r?j), the time 
for the number of heavy stars within the initial half-mass 
radius to increase by 1.5. Here again, the trend is clear: 
the closer the mass ratio is to unity, the longer it takes the 
heavy tracers to sink into the core. 

The shorter core collapse time scales for the models with 
larger /x or larger Wq give these clusters little time to evolve 
dynamically. Consequently, mass segregation is not as pro- 
nounced in these clusters, as is evident from Table 2. 



5. SUMMARY AND DISCUSSION 

We have studied the mass segregation of both light and 
heavy tracers in two-component star clusters. The calcula- 
tions were performed using our recently developed Monte 
Carlo code using N = 10^ stars, with a few key calcu- 
lations repeated using a direct iV-body code running on 
GRAPE computers. The two methods showed good agree- 
ment in all cases. The large that can be used in Monte 
Carlo simulations is essential for this type of study, as 
it allows us to treat a small tracer population of objects 
without suffering large numerical noise. 

We found that light tracers with ^ < 0.1 — 0.4 (de- 
pending on cluster parameters) are ejected from the clus- 
ter core, on average, within one central relaxation time 
(~ lO'' — 10^ yr for most globular clusters). The cores 
of globular clusters should therefore be largely devoid of 
stars of mass < 0.25 Mq. As a result, we predict that any 
observed low-mass star in a dense globular cluster core is 
most probably part of a binary system with a more massive 
companion. For example, low-mass helium white dwarfs 
observed in a cluster core are most probably in binary sys- 
tems (Taylor et al. 2001). 

Some of the low-mass objects, when they are ejected 
from the central region of the cluster, settle in the outer 
halo, where the relaxation time is so long that they are 
prevented from further segregation or evaporation through 
the tidal boundary. We found in our simulations that the 
number density of light tracers in the outer parts of the 
cluster can actually become significantly higher than it was 
initially. We studied the implications of these results for 
gravitational microlensing in globular clusters and found 
that low-mass objects in the cluster halo can dominate the 
optical depth for microlensing by up to an order of magni- 
tude. We applied these results to the recent null detection 
of Sahu, Anderson, & King (2001) to estimate an upper 
limit on the current mass fraction of M22 in very low-mass 
objects of ~ 25% at the 63% confidence level. 

For clusters with heavy tracers (/z > 1) we found good 
agreement betwecin our Monte Carlo results and direct 
A^-body integrations, as well as simple theoretical esti- 
mates. In particular, we found that the time scale for mass 
segregation, tq, varies as l//i, as expected from theoretical 
predictions. Specifically, over a wide range of initial clus- 
ter concentrations and mass ratios, To/tr = C/fi, where 
the constant C ~ 0.5 — 1. 
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